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1. Introduction 

In previous papers [1,2], we extracted the light-cone wavefunctions of the p meson and 
the corresponding twist-2 Distribution Amplitude (DA) using the precise HERA data [3,4] 
on diffractive p photoproduction. In reference [1], we found that the extracted twist- 
2 DA, which is sensitive to the longitudinal light-cone wavefunction, is broader than the 
asymptotic form z(l — z) and agrees very well with QCD Sum Rules predictions. We showed 
in reference [2] that this conclusion is not much model dependent. On the other hand, we 
also found that the data prefer a transverse wavefunction with end-point enhancement [1] 
although the degree of such an enhancement is model dependent [2] . 

In this paper, we extend our study in two ways. First, we shall repeat our fits using 
a CGC-inspired dipole model with a more realistic description of saturation effects than 
the dipole models we used previously [1,2]. Second, we extend the analysis beyond twist-2 
by extracting also the twist-3 DA which is sensitive to the scalar part of the transverse 
light-cone wavefunction. We shall then compare our extracted twist-2 and 3 DAs with 
lattice [5] and QCD Sum Rules predictions [6]. 

In references [1,2], we fit the data using models for the dipole cross section which 
have been previously constrained by F2 data and which also give a good description of 
the diffractive structure function (F®^) data [7,8]. Such models may be called "forward 
models" since, by the optical theorem, the dipole cross section is equal to the forward 
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dipole-proton elastic scattering amplitude. Such models do account for saturation, but in 
a way independent of the momentum transfer t or the impact parameter b. 

More realistic dipole models take into account the t or 6-dependence of saturation 
effects and have been proposed in references [9-11]. Here, we shall use the model by 
Marquet, Peschanski and Soyez [10] and refer to it as the t-CGC model. Given the limited 
Q 2 € [0, 41] GeV 2 range of the HERA data on diffractive p production to which we fit, we 
do not consider DGLAP evolution [11,12]. 

2. Amplitudes and cross-sections 

In the dipole model [13,14], the imaginary part of the amplitude for diffractive p production 
is given by [10]: 

W A ( S , t; Q 2 ) = E / d 2 r dz ^(r, z; Q 2 )*g(r, z)*e~ izr ' A Af (x, r, A) , (2.1) 

h,h 

where t = — | A| 2 . ^ j^( r , z \ Q 2 ) an d Vff'r (r, z) denote the light-cone wavefunctions of the 
photon and p-meson respectively while Af(x, r, A) is the dipole-proton scattering ampli- 
tude. 

We choose the dimensionless variable x carrying the energy dependence of the dipole- 
proton elastic amplitude M(x, r, A) to be a modified Bjorken-x: 

/ 4m 2 A 

X = XB n 1 + -Q^J (2 - 2) 

with a phenomenological quark mass of m/ = 0.14 GeV. This choice coincides with xsj at 
high Q 2 and replaces the latter at low Q 2 . We should point out that the optimal choice 
for the dimensionless variable that carries the energy dependence of the dipole-proton 
amplitude remains an open and interesting problem [15]. For instance, it can be argued 
that the dipole-proton amplitude should not depend at all on the photon's virtuality [16,17]. 

Assuming that J\f(x,r, A) = M(x,r, A), the angular integration can be carried out 
giving 

%mA x (s,t;Q 2 ) = 2ir^2 J dr dz rJ (rzA) f^V, z; Q 2 )^(r, z)*M(x, r, A) . (2.3) 

h,h 

Note that by setting t = in equation (2.1), we recover the forward amplitude used in 
references [1,2], i.e. 

QmA x (s,t;Q 2 )\ t=0 = sY^ [ d 2 r dz *g(r, z; Q 2 )a(x, z)* (2.4) 

h,h 

where we have introduced the dipole cross-section 

£( z> r ) = A/ ~ (X,r, ° ) = - / d 2 b M(x, r, b) . (2.5) 



- 2 - 



This means that the dipole cross-section is independent of the impact parameter 1 and that 
any modelling of saturation at the level of this cross-section is necessarily ^-independent. 
Nevertheless, the dipole cross-section can be extracted using the high quality F2 data and 
can then be used to predict the imaginary part of the forward amplitude for diffractive p 
production and thus the forward differential cross-section: 



dax 
dt 



= -L(WIaM)) 2 (1 + /3 a (0) 2 ) , (2.6) 
t=o 167r 



where f3\ (0) is the ratio of real to imaginary parts of the forward amplitude and is computed 
as in reference [1]. The t-dependence can be restored as suggested by experiment [3]: 



dcf\ _ d<j\ 
dt dt 



x exp(-5|t|) (2.7) 

=0 



where 



1 GeV 2 



0.2 



B = N [ U -°{WTM]) +1 ) (28) 

with N = 0.55 GeV~ 2 . Note that in doing this, we are assuming that it is a reasonable 
approximation to consider the ratio of real to imaginary part of the the amplitude to be 
independent of t, i.e. Px(0) « /3\(t). 
In the non-forward case, 

^ = ±-^ m A x (s,t)) 2 (l + /3 A (t) 2 ) , (2.9) 

where f3\(t) is now the ratio of real to imaginary parts of the non- forward amplitude 
and is computed in the same way as [3\ (0). This differential cross-section is measured at 
HERA [3,4]. 

The total cross-section, which is also measured at HERA, is obtained by integrating ei- 
ther (2.7) or equation (2.9) over t. In the former case, this integration is trivially performed 
analytically for each polarisation and the total cross-section is taken to be a = ol + eaj-? 
In the non-forward case, we integrate equation (2.9) numerically over t to obtain the total 
cross-section. 



3. Dipole models 

In principle, the dipole-proton scattering amplitude M(x, r, b) can be obtained by solving 
the Balitsky-Kovchegov (BK) equation [18-20] which itself can be derived within the Colour 
Glass Condensate (CGC) formalism [21-25]. However, work is still in progress to implement 
in a satisfactory way the impact-parameter dependence in the BK equation [26]. On the 

The dipole-proton amplitude AT and dipole cross-section a are sometimes used synonymously in the 
literature and one speaks of the fe-dependence of the dipole cross-section [11]. 
2 To compare with the HERA data, we take e = 0.98. 
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other hand, an approximate solution for the forward dipole-proton amplitude has been 
proposed several years ago by Iancu, Itakura and Munier [27]: 

x oL | ln(2/rQ s ) 1 
( TO \ 2 | 7s+ KAln(l/a:) I 

Af(rQ s ,x,0) = AU for rQ s < 2 

= {1 - exp[-aln 2 (6rQ s )]} for rQ s > 2 , (3.1) 

where the saturation scale Q s = (xq/x) x ^ 2 GeV. The dipole cross-section is then given 
by <j(rQ s ,x) = aoAf(rQ s ,x,0). In reference [27], the anomalous dimension j s was fixed 
to its LO BFKL value of 0.63 while the free parameters do, xq and A are determined 
by fitting to the F2 data, without a charm contribution. When charm was included in 
reference [11], it was found that the saturation scale decreases dramatically. We used this 
fit in reference [2] where we referred to it as CGC[0.63]. However, Soyez later showed that 
letting the anomalous dimension 7 S vary freely also gives a good F2 fit with no significant 
decrease in the saturation scale [28]. The resulting value of j s is 0.74, close to the value 
obtained with the RG-improved NLO BFKL kernel [28]. Kowalski and Watt confirmed 
the Soyez fit in reference [9] and we used their fitted parameters in reference [2] where we 
referred to it as the CGC[0.74] fit. This model has been extended by Marquet, Peschanski 
and Soyez to include a t-dependence in the saturation scale resulting in what we now refer 
to as the t-CGC model. In the t-CGC model [10], 

Af(rQ s (t),r,t) =jV(rQ s (i),r,0) x exp(-B\t\) (3.2) 

where 

Q s (t) = (x /x)^ 2 (l + c^/\i\) . (3.3) 

The t-CGC model has two additional free parameters (c and B) compared to its forward 
counterpart, i.e. CGC[0.74]. In reference [10], these parameters were fixed by a global fit 
to a subset of the vector meson (p,(j> and J /^) production data including the earlier p data 
from HI [29]. The authors of reference [10] selected those data that were were not too 
sensitive to the precise shape of the meson wavefunction and did not include the CL/cr 
data. As shown in table 1, all other parameters, including the anomalous dimension y s , 
have the same values as in the forward case. This ensures that the quality of the F2 fit is 
unaltered in the t-CGC model compared to that with the CGC[0.74] model. 

The CGC models are based on perturbative QCD and, for that reason, the original 
fits to i*2 were in the range Q 2 G [0.25,45] GeV 2 , thereby excluding the available data 
in the low Q range, Q 2 G [0.045,0.20] GeV 2 , a kinematic region where perturbation 
theory is not reliable [9]. In these fits the light quark mass was neglected compared to 
Q 2 so that the dipole-proton amplitude is evaluated at x = XBj when computing the light 
quark contribution to the structure function [9]. By using equation (2.2), we shall dare to 
extrapolate the use of the CGC dipole models to low Q 2 , including Q 2 = 0. In order to 
legitimate such a procedure, we verify that these models still describe the F2 data in the 
extended range Q 2 G [0.045,45] GeV 2 , where the numerical difference between x given by 
(2.2) and xbj can be significant. The resulting x 2 /d.p values 3 thus obtained are shown in 

3 Hereafter x 2 /d.p stands for xVdata points. 
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the last column of table 1. They indicate that it is reasonable to extrapolate the use of the 
CGC models to low Q 2 , together with the modified Bjorken-x given by equation (2.2). 

Finally, we shall also use the Regge-inspired FSSat dipole model [30], for which the 
original F 2 fits were inclusive of the low Q 2 data, i.e. with Q 2 G [0.045,45] GeV 2 . This 
model falls in the category of forward models, like CGC[0.63] and CGC[0.74], where satu- 
ration effects are assumed to be independent of momemtum transfer. 

Note that we are using here simple phenomenological parametrizations for the dipole- 
proton scattering amplitude subject to the requirement that they fit the F 2 data. Such 
simple parametrizations are convenient for extracting the meson light-cone wavefunction. 
We do not consider here the more sophisticated solutions to the BK equation with running 
coupling [31,32] or the dipole cascade model [33]. 



CGC F 2 fits 



Model 


Is 


x 


A 


<7 [mb] 


B GeV" 2 


c GeV- 1 


X 2 /133 d.p 


X 2 /156 d.p 


CGC[0.63] 


0.63 (fixed) 


1.09 x 10" 7 


0.170 


37.4 






119/133 


135/156 


CGC[0.74] 


0.74 


1.63 x 10" 5 


0.216 


27.4 






113/133 


145/156 


t-CGC 


0.74 


1.63 x 10" 5 


0.216 


27.4 


3.754 


4.077 


113/133 


145/156 



Table 1: The fitted parameters of the CGC dipole models obtained from fits to F2 with Q 2 G 
[0.25,45] GeV 2 and £Bj < 0.01. For each model, the first % 2 /d.p value is obtained when the model 
is fitted to F 2 data with Q 2 G [0.25, 45] GeV 2 . The second x 2 /d.p value is a prediction of the fitted 
model when tested in the extended range Q 2 G [0.045,45] GeV 2 . 

All the forward models we consider also give a good description of the diffractive 
structure function (F®^) data [7,8]. The non forward models have been applied to describe 
a number of processes [34,35] including (in the approach of reference [9]) inclusive hadron 
production, for which the first LHC data are now available [36,37]. 



4. Light-cone wavefunctions 

By analogy with the photon's light-cone wavefunction, the p meson wavefunction can be 
factorised into a spinor and scalar part [38]. In momentum space, 

^i(^ k ) = \/§^^ k )^(^k) (4.1) 



where 



a,.. %(^+-k) x v h ((l-z)P + ,k) 



S p 'Uz,k)= ' V .7 , ^ ' (4.2) 

where P + is the "plus" component of the 4-momentum of the meson given by 

p " = r + '27^' 0± ) • (4 - 3) 
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Using the longitudinal polarisation vector 



M~ p ' 2P+' ± 



we can write 



where [39] 



and 



so that 



Mp 
2P^ 



-M p 8 h - h 



PS* (*>*)]- 



p+ 



- v _h_ 



mj + k 2 

,2(1 - z)M n 



J h,-h 



M p z{l - z) 



[z(l-z)M 2 p +m} + k 2 ] 5 h; _ n . 



To evaluate the transverse spinor wavefunction, we use the polarisation vector 

1 



to obtain [39] 



V2 



V2 



(0,0, e±) 



with 



(l,±i) 



(4.4) 

(4.5) 
(4.6) 

(4.7) 
(4.8) 

(4.9) 



{[^h±,h T " (1 " z)]6^- h± ke ±ie « + m/^si} (4.10) 



h.h 



z(l-z) 



where we have used the polar representation of the transverse momentum, i.e. k = ke lSk . 

Taking the two-dimensional Fourier transforms of equation (4.1), we obtain the r-space 
light-cone wavefunctions [38]: 



*) = -V 2 6h '-' h M pZ {i-z) [z{l - z)M ? + m f~ v ^ i(r ' z) (411) 




1 



where V 2 = ^d r + d 2 , and 



P.T(7=±). 
h,h 



r,z) = ± 




47T z(l - z) 



[ie 1 r (z5 h± - hT - (1 - z)5 hTi ft ± )<9r + mf5 h ±fi ± ](j)T(r,z). 

(4.12) 

Previous work [9,10,33,38] has argued that a reasonable assumption for the scalar part of 
the light-cone wavefunction for the p is of the form 



'm)R\ s 



m)R\ 



<^ G (r,z) = N\ 4[*(1 - z)] b y2vri?2 exp 1^1 exp I - g [jg(1 V _ ^ I 



x exp 



2[z{l - z)] 6 *r 2 
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which is referred to as the "Boosted Gaussian" (BG) wavefunction. 4 These light-cone 
wavefunctions are normalised: 

J d 2 rdz\V p > x (r,z)\ 2 = 1 , (4.14) 

where we have defined 

\^(r,z)\ a = ^\9^r,z)\ 2 . (4.15) 

h,h 

This normalisation condition embodies the assumption that the p consists of the leading 
qq Fock state only. In addition, the longitudinal light-cone wavefunction is subject to the 
leptonic decay width constraint [38]: 



V2f p M p = ^ I dz [z(l - z)M 2 p + m) - V 2 ]% 

TT J Z(l 



4>L(r,z) 



z) 



(4.16) 

r=0 



where the decay constant f p is simply related to the experimentally measured electronic 
decay width T p _^. e + e - of the p [1]. 

5. Fitting the HERA data 

In the original BG wavefunction, bi = &t = 1 and B? L = = 12.9 GeV -2 so that 
the leptonic decay width and the normalisation constraints are satisfied. However, when 
this BG wavefunction is used in conjunction with either the FSSat or any of the CGC 
models, none of them is able to give a good quantitative agreement with the current 
HERA data. This situation is considerably improved by letting R\ and b\ vary freely and 
this enhances end-point contributions. We also investigate the requirement for additional 
end-point enhancement in the transverse wavefunction by using a scalar wavefunction of 
the form 

Mr, z) = 4>T G (r, z) x [1 + c T e + d T ^} (5.1) 

where £ = 2z — 1. 

We fit to the same data set as in reference [1], i.e. to total cross-section data, the ratio 
of longitudinal to transverse cross-section data and the decay constant datum for the longi- 
tudinally polarised meson, i.e. to a total of 76 data points. For the non-forward models, we 
also include the differential cross-section data with \t\ < 0.5 GeV 2 (46 data points) result- 
ing in a total of 122 data points. As in reference [1], we rescale the overall normalisation 
of the data down by 5% which is consistent with the experimental uncertainty. 

Note that for the non-forward t-CGC model, we make use of the current data on the 
differential cross-section to refit the parameters controlling the t-dependence. This means 
that there are 2 additional free parameters, B and c, compared to the forward FSSat, 
CGC[0.63] and CGC[0.74] fits. 

Our best fit parameters and resulting x 2 /d.o.f are given in table 2. For the t-CGC 
model, we find two successful fits which we refer to as t-CGC and t-CGC(alt.). We choose 



4 This wavefunction is a simplified version of that proposed originally by Nemchik, Nikolaev, Predazzi 
and Zakharov [40]. 
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to record the alternative t-CGC fit because we extract from it a qualitatively different 
shape for the transverse wavefunction: the parameters cx and dx go negative in this fit. 

We note that our fitted values for the parameters controlling the t-dependence are 
B = 3.159(3.229) GeV" 2 and c = 1.146(1.262) GeV" 2 , where the number in brackets 
corresponds to the alternative t-CGC fit. The values for c are lower than that obtained in 
reference [10] and which we quote in table 1. 

Our fits to the HERA data are shown in figures 3 and 4 for the total cross-section, in 
figure 5 for the gl/ot ratio. In figures 6 and 7 we show the t-CGC fit to the differential 
cross-section data as well as the predictions of the forward models FSSat, CGC[0.63] and 
CGC[0.74]. 



Best fits for the meson wavefunction 



Reference 


Model 


R\ 


R^ 


b L 


bx 


c T 


dx 


X 2 /d.o.f 


[1,2] 


FSSat 


26.76 


27.52 


0.5665 


0.7468 


0.3317 


1.310 


68/70 


[2] 


CGC[0.63] 


27.31 


31.92 


0.5522 


0.7289 


1.6927 


2.1457 


67/70 


[2] 


CGC[0.74] 


26.67 


21.30 


0.5697 


0.7929 








64/72 


This paper 


t-CGC 


29.63 


21.64 


0.5018 


0.7396 








114/116(72/75) 


This paper 


t-CGC (alt.) 


29.68 


20.96 


0.5004 


0.7314 


-0.1601 


-0.1656 


112/114(70/75) 



Table 2: Best fit wavefunction parameters and resulting x 2 /d.o.f obtained with each dipole model. 
For the non-forward t-CGC model, the partial x 2 /d.p for the total cross-section and ratio data is 
given in brackets. 

In figures 1 and 2, we show the extracted light-cone wavefunctions squared, as defined 
by equation (4.15), and evaluated at r = 0. Figure 1 confirms that the longitudinal 
wavefunction is not much model-dependent. On the other hand, figure 2 shows that the 
data allow for a wide range of end-point enhancement in the transverse wavefunction. 

Additional theoretical constraints on the p light-cone wavefunction can be obtained 
from QCD Sum Rules and lattice QCD. More precisely, these non-perturbative methods 
can predict the moments of the DA at a low non-perturbative scale. In order to use these 
constraints we must first establish the relationship between the DA and the light-cone 
wavefunction. We shall then be in a position to compare the DA predicted by each of our 
extracted wavefunctions with QCD Sum Rules and lattice predictions. 



6. Distribution Amplitudes 

Light-cone DAs appear in the operator product expansion of vacuum-to-meson transition 
matrix elements of quark- antiquark non-local gauge invariant operators at light-like sepa- 
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0.2 0.4 0.6 0.8 1 

z 

Figure 1: The longitudinal light-cone wavefunction squared evaluated at r = 0. Solid: CGC[0.74], 
Dotted: FSSat, Dashed: CGC[0.63], Dot-dashed: t-CGC. The alternative t-CGC(alt.) fit is indis- 
tinguishable from the t-CGC fit. 

0.03 



0.025 



0.02 

-S 
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"0 0.2 0.4 0.6 0.8 1 

Z 

Figure 2: The transverse light-cone wavefunction squared evaluated at r = 0. Solid: CGC[0.74], 
Dotted: FSSat, Dashed: CGC[0.63], Dot-dashed: t-CGC, dash-dash-dotted: t-CGC(alt.). 

rations: 5 [44] 

(0|g(0)[0,x] 7 ^(x)|p(P,A)) = V2f p M p e -^P» C &u e~ iuP ' x 4>\\ (u, fi) 

" " x Jo 

+ (fi ~ P^jr^) ^fpM P f du e- luP - x g ± (u, M ) 

+ {•••} ° (6-1) 

5 We use the classification and notation of Ball and Braun [41, 42] with the notational simplication 
g\ — > g± since we consider only a vector coupling here. In the alternative notation of reference [43], 
4>\\ = <pi and g± = ip^. We also note that in references [41,42] the decay constant is denned as y2f p . 
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Figure 3: 

CGC[0.63], 
since at hig 



Fits to the HERA cross-section data. Solid: CGC[0.74], dotted: FSSat, dashed: 
dot-dashed: t-CGC. The alternative t-CGC(alt.) fit is only shown at Q 2 
her Q 2 it is indistinguishable from the t-CGC fit. 



GeV^ 



3.2) 



where [i is the renormalization scale, x ■ x = 0. The gauge link 

[0,x] = Pexp[ig[ dt(-x)„A»({l-t)x)], 
Jo 

ensures the gauge invariance of the matrix elements. According to the classification of Ball 
and Braun [42,44], <f>» is twist- 2 and g± is twist-3. The {• • • } in equation (6.1) stand for 
higher twist contributions [42,44] which we do not consider here. 
The //-dependent DAs 

<p = {<f>\\,g±} (6.3) 

satisfy the following normalisation condition: 

-l 

du <p(u, ju) = 1 . (6-4) 



I 

Jo 



This guarantees that in the local limit x — > 0, (6.1) becomes 

<0|?(0) 7 "g(0)|p(P,A)> = V2f p M p e» 



(6.5) 
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80 120 

W [GeV] 



Figure 4: Fits to the ZEUS total cross-section. Solid: CGC[0.74], dotted: FSSat, dashed: 
CGC[0.63], dot-dashed: t-CGC. The alternative t-CGC fit is only shown at Q 2 = 0.47 GeV 2 
since at higher Q 2 it is indistinguishable from the t-CGC fit. 



thereby recovering the definition of the decay constant f p . 

In order to establish a connection between these DAs and the light-cone wavefunctions 
of the p, we need to apply the equal light-cone time condition, x + = and choose the 
light-cone gauge A + = 0. 6 We can therefore rewrite equation (6.1) as 

(0|g(0) 7 ^(0,x-,0 ± )|p(P,A)) = ^/ P M^P" / tin e~ iu p+x ~ ^(u, p) 

+ V2f p M p (< - P"^^) £du e- iuP+x -g ± ( U ,p)(6.6) 

To isolate <f>», we take \x = + with A = L in equation (6.6) and obtain 

(0\q(0)l + q(^)\ P (P,L)) = V2f p P+ J due-^-^u^) . (6.7) 

6 It is under these conditions that the Fock expansion of the physical p is carried out [39]. 
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Figure 5: The a L /cr T data are at W = 75 GeV for HI and W — 90 GeV for ZEUS. 
Solid: CGC[0.74], dotted: FSSat, dashed: CGC[0.63], dot-dashed: t-CGC, dash-dash-dotted: t- 
CGC(alt.). 



On the other hand, to isolate g±, we take the scalar product of equation (6.6) with the 
complex conjugate of the meson's transverse polarization vector to arrive at 

(0|?(0)4 • yq(x-)\p(P,T)) = V2f p M p j du e~ iuP+x ~ g ± (u, M ) . (6.8) 

Taking the Fourier transform of the above matrix elements with respect to the longitudinal 
light-cone distance x~ yields: 

<j>\\(z,n) = ^rj dx-e* zP+ *- (0\q(0h + Q(x-)\p(P,L)) (6.9) 

and 

g L (z,f,) = [ dx-e izP+ *-{0\q(0)e T *.7q(x-)\ P (P,T)) . (6.10) 

y/2jpMp J 
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Figure 6: Dot-dashed: t-CGC fit to the differential cross-section data from ZEUS in the range 
< |t| < 0.5 GeV 2 . At these values of Q 2 the alternative t-CGC(alt.) fit is indistinguishable from 
the t-CGC fit and is not shown here. Also shown are the predictions of the forward models: Solid: 
CGC[0.74], dotted: FSSat, dashed: CGC[0.63]. 



We now wish to express the matrix elements on the right-hand-sides of the above 
equations (6.9) and (6.10) in terms of the meson's light-cone wavefunction given by equation 
(4.1). To do so, we use the Fock expansion 



h.h 



dfc+d 2 k 
16ttV^ + (^ + - k + ) 



^ h (k + /P + , k)P h (k+, k)4((p+-* : +) > -k)|o) 



(6.11) 



and the mode expansion 

, "dA; + d 2 k 

g(0,x-,0) 



J2K(k+, k)u a (k+, k)e~ ik+x ~ + di(k+, k)v a (k+,k)e 



ifc+a 



167T 3 /C+ 

together with the equal light-cone time anticommutation relations 

{b[(k + ,k),b h/ (k + ',k')} x+=0 = 167r 3 k + 5(k + - k + ')5 2 {k-k')5 hh , etc, 



(6.12) 
(6.13) 
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Figure 7: Dot-dashed: t-CGC fit to the differential cross-section data from HI in the range 
< \t\ < 0.5 GeV 2 . At these values of Q 2 the alternative t-CGC(alt.) fit is indistinguishable from 
the t-CGC fit and is not shown here. Also shown are the predictions of the forward models: Solid: 
CGC[0.74], dotted: FSSat, dashed: CGC[0.63]. 



in order to express the meson-to-vacuum matrix elements as 



(0|S(0)yg(aT)|p(P,A)) = y^N c Y, I 

h,h 



dfc+d 2 ke(|k| < fi) 



¥$(fc+/p+k) 



x v- h (P + -k + ,-k)^u h (k + ,k)e 



ik + 3 



(6.14) 



Here we have identified the renormalization scale /i with the ultraviolet cut-off on transverse 
momenta [45,46]. 

Taking the Fourier transform with respect to the light-cone distance x~ allows us to 
carry out the integral over k + , which in turn fixes the plus- momentum of the quark to be 
zP + and so 



P + I <L ! ;"e-" 2P+ (0|g(0) 7 ^(x-)|p(P,A)) = ^ ^*g(*,k) ( 6 - 15 ) 



^((l-^P+.-k) „u fc (zP+,k) 



We can now take the scalar product with the meson's transverse polarisation vector to 
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obtain 

/ p\k\<fi, j2i_ 

dx-e™- zP+ (0\q(0)e T *.jq(x-)\ P (P,T)) = N c ^ J k)| 2 r (z, k) 

h,h 

(6.16) 

and with the help of equation (6.10), we can deduce that 

^ 2 -" ) " TlfiS E S'^'^^.k) . (6.17) 

Jp p h,h 

On the other hand, if we consider A = L and n = +, equation (6.15) becomes 

//■|k|</i ^21, 
dx-e^-^^l^+ff^-JIpCP.L)) = -2iV c ^ / — 3 S^ h (z,k)Mz,V5 h ^ 

h,h 

(6.18) 

and so 

^ z, ^ = ~V^f~p^J T^ s hi ( - z,k ^ z,k)6h ^-~ h ' (6 ' 19) 

Using equations (4.8) and (4.10), we can rewrite (6.17) and (6.19) as 

^^ ) = wf^y w [m / + (z +(1 - z))k] ^r^F (6 - 20) 

and 

* "> = ^ (W [M > (1 " + *"? + k2] (6 ' 21) 

respectively. 

Finally, inserting the Fourier transform pairs 

j d 2 re- !k ^( Z ,r) = ^(z,k) (6.22) 

and 

y d 2 re-* k - r {-V 2 A (z, r)} = k 2 A (z, k) (6.23) 
in (6.20) and (6.21) and carrying out the integration over |k|, we arrive at 

(*,/.)= *!* [ drM(H[M 2 z(l-z)+m 2 -V 2 ]^4 ^ 

and 

" sfe / dr " Jl( " r) [(m? " (z2 + (1 " 2)2,v ' ] ^? ' (6 ' 25) 

Equations (6.24) and (6.25) are our main results. They show that the DAs can be expressed 
in terms of the scalar parts of the light-cone wavefunctions. 
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Note that in the limit \x — > oo, equation (6.24) becomes 

N 



\(z,oo) 



(6.26) 



r=0 



so that the decay width constraint given by equation (4.16) implies that our extracted twist- 
2 DA 0|i satisfies the normalisation condition given by equation (6.4). On the other hand, 
our extracted twist-3 DA g± satisfies only approximately this normalisation condition. 



7. Comparison with Sum Rules and lattice predictions 

To make contact with Sum Rules predictions, we expand the twist-2 DA in Gegenbauer 
polynomials, i.e. [44] 



<t>\\(z,n) = 6z(l - z) 



i=2,4,... 



(7.1) 



where C- are the Gegenbauer polynomials. The factor z(l — z) is the weight function 
over which the Gegenbauer polynomials are orthogonal. Keeping only the lowest conformal 
spin contribution in the expansion yields 



<t>\\(z,n) = 6z(l - z) 



(7.2) 



Similarly an explicit expression for the twist-3 DA is [44] 



+ 



il2 a 2^) + 64 ( 3w 3 , (^) ~ 



(3 - 30^ 2 + 35^ 4 ) 



(7.3) 



The QCD sum rules estimates for the various non-perturbative parameters have been 
updated in reference [6]: a\ = 0.15 ± 0.07, C3 = 0.030 ± 0.010, oj% = 5.0 ± 2.4 and 
ojf = —3.0 ± 1.4, at a scale of \i = 1 GeV. Their perturbative evolution with the scale /j, 
are explicitly given in reference [44]. As \x — > 00, all these parameters vanish and we obtain 
the asymptotic DAs given by 

<j>\\(z,oo) = 6z(l-z) (7.4) 

and 



g±(z, 00) = 7 (l + £ 2 



(7.5) 



We are now able to compare our extracted DAs computed using equation (6.24) and 
(6.25) with the QCD Sum Rule prediction given by equation (7.2) and (7.3) respectively. 
We first confirm our earlier observation [1] that our extracted DA hardly evolves with the 
scale fx when fj, > 1 GeV, i.e. it neglects the perturbatively known /x-dependence. Given 
the limited Q 2 range of the HERA data to which we fit (y^Q 2 < 7 GeV), our extracted 
DA should thus be viewed as a parametrization at some low scale /i ~ 1 GeV. 
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Figure 8: The extracted twist-3 DAs compared to the QCD Sum Rules predictions. Dotted: FSSat; 
Solid: CGC[0.74]; Dashed: CGC [0.63]; Dot-dashed: t-CGC; Dash-dash-dotted: t-CGC (alt.); Long- 
dashed: Sum Rules at 1 GeV; Dot-dot-dashed: Sum Rules at 3 GeV. Crosses: Asymptotic. 



We compare in figure 8 our extracted twist-3 DAs with the Sum Rule DA at two values 
of ji: 1 and 3 GeV. Note that we obtain the DA at 3 GeV by evolving the DA at 1 GeV 
to leading logarithmic accuracy [44]. All the distributions show an enhanced end-point 
compared to the asymptotic DA given by equation (7.5) and which is also shown on the 
plot. Note that the non-monotonic behaviour of the Sum Rule distribution is not physical 
and is due to the truncation in the Gegenbauer expansion. 

In figure 9, we show our extracted twist-2 DAs which are sensitive to the longitudinal 
wavefunction, i.e. computed using equation (6.24). Similar predictions appeared in ref- 
erence [1,2] for the FSSat and the forward CGC models. We now also show the twist-2 
DA extracted with the non-forward t-CGC model and in addition we use equation (6.24) 
to compute the DA instead of assuming that the DA is simply proportional to the scalar 
light-cone wavefunction integrated over transverse momentum [47]. 

The extracted DAs are all consistent with the QCD Sum Rule prediction at 1 GeV, 
given the size of the uncertainty on the parameter a\ ■ Table 3 compares the corresponding 
second moments of the twist-2 DA given by 

(e 2 ) M = / d*£Vu (*,/*). (7.6) 

J 

We also show in this table the prediction of the old Boosted Gaussian wavefunction which 
does not fit the HERA data. 

8. Conclusions 

We have extracted the leading twist-2 and subleading twist-3 Distribution Amplitudes of 
the p meson using the HERA data. We find that the twist-2 DA is not much model 
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Moments of the leading twist DA at the scale fj, 



Reference 


Approach 


Scale n 




[1] 


Old BG prediction 


~ 1 GeV 


0.181 


This paper 


FSSat fit 


~ 1 GeV 


0.267 


This paper 


CGC[0.74] 


~ 1 GeV 


0.266 


This paper 


CGC[0.63] 


~ 1 GeV 


0.271 


This paper 


t-CGC 


~ 1 GeV 


0.286 


[5] 


Lattice 


2 GeV 


0.24(4) 


[6] 


Sum Rules 


1 GeV 


0.254 


[6] 


Sum Rules 


3 GeV 


0.237 




6z(l - z) 


oo 


0.2 



Table 3: Our extracted values for (^ 2 ) M , compared to predictions based on the QCD Sum Rides 
or lattice QCD. Note that the Sum Rules prediction at fj, = 1 GeV has an uncertainty of ±0.024 
corresponding to the uncertainty of 0.07 in a\. 
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Figure 9: The extracted leading twist- 2 DAs compared to the Sum Rule predictions. Dotted: 
FSSat; Solid: CGC[0.74]; Dashed: CGC[0.63]; Dot-dashed: t-CGC; Long-dashed: Sum Rules at 1 
GeV; Dot-dot-dashed: Sum Rules at 3 GeV. Crosses: Asymptotic. 



dependent where as the data allow for a family of extracted twist-3 DAs with a varying 
degree of end-point enhancement. The extracted DAs are consistent with the QCD Sum 
Rules and lattice predictions. 

In its present form, our parametrization for the DA lacks the perturbative QCD evo- 
lution. Improving this would enable a more precise comparison with Sum Rules and the 
lattice. 
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